Moving glass phase of driven lattices 
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We study periodic lattices, such as vortex lattices, driven by an external force in a random 
pinning potential. We show that effects of static disorder persist even at large velocity. It results in 
a novel moving glass state which has analogies with the static Bragg glass. The lattice flows through 
well-defined, elastically coupled, static channels. We predict barriers to transverse motion resulting 
in finite transverse critical current. Experimental tests of the theory are proposed. 
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' An open question is to understand the effect of static 
, substrate disorder on periodic media such as vortex lat- 
CNj ■ tices p|, charge density waves (CDW) B], Wigner crys- 
tals [|[^ colloids Q, magnetic bubbles Numerous ex- 
' periments study such elastic systems in motion under an 
, applied force produced by a current (vortex lattices), a 

■ voltage (CDW), an electric field (colloids) and a mag- 
i netic field gradient (magnetic bubbles). It is therefore 

■ important to describe the physical properties of both the 
static and moving lattices. The statics of vortex lattices 
has been much investigated recently and it is generally 
agreed that disorder leads to a glass state with diverg- 
ed ' ing barriers, pinning and loss of translational order. The 
^ ] precise nature of the glass state however has been the 

I ■ subject of much debate in particular concerning the 
' ^ , decay of translational order and the presence of topolog- 
^ ' ical defects. It was shown recently ||] within an elastic 
Q ] theory that, because of the periodicity of the lattice, the 
decay of translational order is only algebraic and that 
. !^ I the resulting glass phase still exhibits divergent Bragg 
' peaks. We argued that at weak disorder, a Bragg glass 
^ ] exists without equilibrium dislocations and also that such 
- - ' a glass will undergo a transition into a strongly disordered 
vortex glass containing topological defects, or a pinned 
liquid, upon increase of disorder or field This is com- 
patible with recent decoration and neutron experiments 
pO| , pl| and with the behaviour of the critical point ||l^] 
in the phase diagram of vortex lattices where a transi- 
tion between two different glass states is observed upon 
raising the field [ p^ . 

It is thus crucial to determine how much of the glassy 
properties of the static system remain once the lattice is 
set in motion, and how translational and topological or- 
der behave. At large velocity v it was expected that since 
the pinning force on a given vortex varies rapidly, disor- 
der would produce little effect. Perturbation theories in 
disorder and 1/v were thus developped to compute 

velocity as a function of the external force, and to esti- 
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mate critical currents. Recently, Koshelev and Vinokur 
p6t have extended the perturbation theory of to 
compute vortex displacements u. They concluded that 
at low T and above a certain velocity the moving lattice 
is a crystal at an effective temperature T' — T -\- Tgh- 
Several experiments indeed suggest that a fast moving 
lattice is more ordered [|l^,|o| . The effect of pinning can 
be described by some effective shaking temperature 
Tsh l/w^ defined by the relation {\u{q)\'^) = Tsh/ceeq'^ 
This would suggest bounded displacements in d > 2 and 
absence of glassy properties in the moving solid. 

In this Letter we reconsider this problem. We show 
that in the case of a moving lattice the perturbation the- 
ory of breaks down, even at large v. The physical 
reason is that some modes of the disorder are not af- 
fected by the motion and static disorder is still present 
in the moving system. As a result the moving lattice is 
in fact a moving glass. Since translational order in the 
moving frame decays and relative displacements are not 
bounded, such a phase cannot be described by simple 
perturbation theory [p^ . As in the statics, periodicity 
is crucial and the moving lattice has a completely differ- 
ent behaviour than other driven systems such as mani- 
folds. The physics of this new phase can be described 
in terms of elastic channels. When the force is applied 
along a principal lattice direction the rows of the lattice 
flow along well-defined, nearly parallel, preferred paths 
in the pinning potential. The manifold of these optimal 
channels (lines for 2D lattice and sheets for 3D vortex 
line lattice) which exhibits a roughness that we estimate, 
are a purely static and reproducible feature of the disor- 
der configuration. We also predict that the moving glass 
exhibits barriers to an additional small transverse force 
and compute the associated transverse critical current. 
The other modes of the disorder are suppressed by mo- 
tion and give rise to an additional wiggling motion of the 
particles around the static channel configuration, which 
can be treated in perturbation. 

We now derive the equation of motion for a lattice 
submitted to external force F. We denote by Ri{t) the 
position of an individual vortex in the laboratory frame. 
The lattice as a whole moves with a velocity v. We thus 
introduce the displacements Ri{t) — +vt+Ui{t) where 
the R^ denote the equilibrium positions in the perfect 
lattice with no disorder. Ui represent the displacements 
in the moving frame. We consider in the following the 
elastic limit in the absence of topological defects, thus 
assuming Uj+i | <^ a where a is the lattice spacing, an 
assumption which may be checked self-consistently. One 
then takes the continuum limit Ui{t) — > u{r,t), where 
u{r,t) is a smoothly varying n-component vector field, 
which components we denote by Ua{r, t). It is convenient 
here to express the displacement field Uq t) in terms 
of the coordinates (r, t) of the laboratory frame. The 
equation of motion in the laboratory frame is then: 
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rjdtUa + -nv ■ VUa = cV^Ma + fP'" + F„ - rjVa + (a (1) 

where 77 is the friction coefhcient, F^™(r,t) = 
—S£/5ua{r,t) is the pinning force, £[u{r,t)] is the 
pinning energy, and the thermal noise satisfies 
Ca(r, t)C0{r',t') = 2TSapS'^{r - r')S{t - t'). For clarity 
we use here an isotropic elastic constant c. The realis- 
tic case, discussed at the end, has the same large dis- 
tance physics. The term rjv ■ Vmq comes from expressing 
the displacement field in the laboratory frame and —rjVa 
is the average friction, v is determined by the condi- 
tion that the average of u is zero. (1) is exact up to 
higher powers of derivatives of u, negligible in the elastic 
limit. The pinning energy can be expressed in terms of 
the vortex density p{r,t) = J2i^i^ — R° — vt — Ui{t)). 
One has £[u{r,t)] — J dfirp{r^t)V{r) where the random 
potential has correlations {V{r)V{r')) ~ A(r — r') of 
range Vf. Since even for smooth displacement fields 
the density is a series of delta peaks, the continuum 
limit for £[u\ should be performed by distinguishing |^] 
the various Fourier components of the density p{r, t) = 
Pq{1 — V ■ u -\- X^K^^o exp(ii4r ■ [r — vt — u(r, t))) where K 
spans the reciprocal lattice and po is the average density. 
Using this decomposition in (|l|) the force due to disorder 
naturally splits into a static and a time-dependent part: 

TjdtUc. + T]v ■ Vu„ = cV^Ua + "*(r, u{r, t)) + F^y'^(r, t, u{r, t j) + - tjv^ + Ca(r, t) (2) 

F^'-'ir,u) = V{r)po ^ ii^^e^^'t'-") - poV„T/(r), F^y'\r,t,u) = V{r)po ^ zif^e^^' 

K.v=0 K.v=^0 

The static part of the random force comes from the 
modes such that K.v = which exist for any direction 
of the velocity commensurate with the lattice. The max- 
imum effect is obtained for v parallel to one principal 
lattice direction, the situation we study now. Reflection 
symmetry then imposes that v and F are aligned along 
direction x, the d—1 transverse directions being denoted 
by y. F^''^' gives the dominant contribution to the lattice 
deformations. In first approximation we drop f^^" and 
solve the remaining static problem (leading to a reference 
ground state at T = 0). F^y" gives additional fluctua- 
tions on top of this ground state, estimated below. The 
static term which comes from the Fourier com- 

ponents k <^ 1/a of the disorder, produces alone only 
bounded displacements for d > 1. Thus, as for the non- 
moving lattice ^ for d > 2, it does not change the large 
scale physics and we drop it. Since F^'**' is now along y, 
and depends only on Uy, — in the ground state. 

The most important terms in (1) thus lead to the fol- 
lowing equation of motion in the laboratory frame which 
involves only the transverse displacements Uy: 

rjdtUy + rjvd^uy = cV^Uj, -f F"'^*(r, Mj,(r, t)) + Ca(r, 

F'''^\x,y,uy) = V{x,y)po ^ KySmKy{uy - y) (3) 

This is now a non-trivial static disordered model and one 
expects a glass phase at low temperature, with pinning of 
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the field u(r, t) into preferred configurations. Thus, the 
moving vortex configurations can be described in terms of 
static channels that are the easiest paths where particle 
follow each other in their motion. Channels in the elastic 
flow regime behave differently than the one introduced 
to describe slow plastic motion between pinned islands 
p9[ . In the topologically ordered moving glass they form 
a manifold of elastically coupled, almost parallel lines or 
sheets (for vortex lines in d = 3) directed along x and 
characterized by transverse wandering Uy. In the labora- 
tory frame they are determined by the static disorder and 
do not fluctuate with time. In the moving frame, since 
each particle is tied to a given channel which is now mov- 
ing, it indeed wiggles and dissipates but remains highly 
correlated with the neighbors. To obtain the roughness of 
the manifold of channels we compute the correlator of rel- 
ative displacements B{x, y) = {[u{x, y) — u{0, 0)]^). A de- 
tailed analysis will be presented elsewhere Q . One de- 
fines two characteristic lengths for decay of translational 
order, and Ry along the longitudinal and transverse 
direction by B{R) ^ . One expects three regimes. 

short scale regime: At very short scales one can ex- 
pand the pinning force to lowest order in u. This 
gives a simple model where pinning is described by a 
random force F^^'^^{x) independent of u whose corre- 
lator is (F'^tat^^)^stat^^/)^ ^ ^gd(^^ _ ^i^j^ ^ ^ 

Po^Ky ^y^f<:- "^^^^ dynamic equivalent of the 

Larkin random force model and B = Brf + {u'^)th'- 

I N /" dq^d'^'^Qy A(l - cos(g^a; + qyy)) 

and (w^)j^ is the thermal displacement. One finds for 
X > c/rjv B{x,y) ^ IS.^^^^ H {cx / rjvy ) where HiQ) = est 

and H{z) ^ 2;'^^"'')/^ at large z. x scales as y^ and the 
displacements are very anisotropic. For x < c/rjv one has 
the usual isotropic result. If one could extrapolate this 
behaviour to larger scales it would result in an algebraic 
decay of translational order in d — 3 {B{x,y) ^ log \ y\) 
and exponential decay in d = 2. However, since the 
Larkin model rests on the expansion in powers of m, it 
is valid only as long as A'maxM 1, where i^max ^ ^/i~f 
is the highest Fourier component of the random poten- 
tial ^. This defines two lengths i?^ and Ry such that 
B{R) ^ r^, below which the Larkin model is valid. At 
large velocity: R^ = (r^wc/A)!/^^-^)^ ^ v{Rlf/c. 
For smaller velocities v < v* ^ c(A/c^ry)^/'^^~'*-' the elas- 
tic term dominates and i?^ ~ Ry ~ Riso where R^^g is 
the static pinning length. These lengths are renormalized 
by temperature and by the dynamical part u'^. Note that 
this Larkin random force corresponds formally to the so- 
called "random mobility term" considered in (2^,^^ and 
by keeping only this term one misses all the physics of 
the moving glass, e.g. the channels and the transverse 
barriers. As for the static case the pinning force in 
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should be treated to all orders in u. Above this length 
scale pinning and metastability appear. 

intermediate regime: At intermediate scales Ry < y < 
Ry and R% < x < R%, the analogous of a random mani- 
fold regime |^,|| exists for which Uy ~ y^. The channels 
are determined by optimization of elasticity {cq^ term), 
dissipation (irjvqx term) and the random potential seen 
independently by each channel in its vicinity. One ex- 
pects many metastable nearly optimal configurations in 
that regime and glassy behaviour. Flory type arguments 
suggest that the scaling properties of this glass are related 
to the static Bragg glass by d — ^ d + 1 and n n — 1. 
The former comes from assuming ^ qy and the latter 
from u ^ Uy. The Flory estimate is then (-^ = 

asymptotic regime: At large distances x > R%,y > Ry, 
in d = 3 the displacements have a slower, logarithmic 
growth. Estimates a la Fukuyama-Lee then give: 

Rl ^ (a\c/A)i/(3-'^), Rl = v{Rlflc (5) 

The moving glass is highly anisotropic since R'^/Ry di- 
verges as u — > cxD. Its upper critical dimension is d = 3, 
instead of d = 4 for the static one. For d > 3 the moving 
system is not a glass but a perfect crystal at weak dis- 
order or large v. For d < 3 weak disorder destroys long 
range order and results in a moving glass. 

As an important consequence of the existence of the 
moving glass, barriers for transverse motion exist once 
the pattern of channels is established. Thus the response 
to an additional small transverse force Fy is very non 
linear with activated behaviour. At T = 0, neglecting 
the dynamic part of the disorder a true transverse criti- 
cal current Jy should exist. This can be seen by adding 
a transverse force in (||). For r]v < rjv* ^ {Riso/'''f)^c 
where Fc is the isotropic critical force, the Larkin do- 
mains remain isotropic and one expects Jy to decrease 
slowly from Jj? ^ to a fraction of the longitudinal critical 
current Jf^^ (since only the Ky modes contribute). For 
rjv > rjv* one expect a much faster decay and a naive es- 
timate for Jy is obtained by balancing the pinning force 
with the transverse force acting on a Larkin domain: 

where A = A/t; is an effective velocity-dependent dis- 
order and Cl is the speed of light. In d = 3 it yields 
exponential decay with v. In practice, rjv* can corre- 
spond to a large driving compared to Fc- The above 
regimes correspond to collective pinning with Ry > a. 
For Ry < a, i.e v < vo — Aa^^'^ /crj:, single channel pin- 
ning leads to different estimates for Jy. Since motion is 
not modified below Jy, dvx/dfy also vanishes below the 
transverse critical force. One can expect non linear ef- 
fects in the flow along x since channel configuration is 
modified when F^ is increased. A simpler example of 
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transverse barriers is a lattice driven in a commensu- 
rate washboard potential V{x,y) = UoCOs{Koy) — F^x. 
There it is clear that even in the moving frame the prob- 
lem is static and that the transverse critical force is 
Fy ~ UqKq. Even a single particule in the 2D poten- 
tial V{x,y) = fcCos{x) + /cCos(j/) — F^x, has finite 



Its velocity can be computed as in |23| and becomes for 
T^O,Vy^ (F2 - flY'^e{Fy - /J, independent of F,. 

To estimate the effects of the time-dependent pinning 
force in (|^) , we split u = u'^ + u'^ into a static u'^ and a 
dynamics u'^ part. A reasonable estimate is: 

pIk'^/^k5{lo - K. 
rpv'^{K^ + q^Y -t- c^q 



d d-. _ V- ptK'/^K5{uj ^ K.v) 



K.v^O 

The dynamic correlations are bounded due to the pres- 
ence of mass term KxV in the denominator, u"^ saturates 
at large distances, even in d = 2 if T = 0. u"* is smaller 
by a factor a/Rx ^ 1 compared to at the length scales 
Rx- u'^ thus represents a small additional wiggling mo- 
tion around the ground state. The massive propagator 
in ^ is very different from a thermal one l/q^ ■ 

Extension to realistic elastic energy, e.g. a triangular 
lattice in d = 2 is straightforward. The static displace- 
ments, within the random force model, are now Ua{q) = 
Fy{l){Ply{'irivqx+CQQq^)~^+P^y{irivqx+Ciiq^)~^) where 
cii and C66 are (dispersionless) bulk and shear mod- 
uli, respectively. Thus the mean square displacement 
B{x,y) is again given, for y > y* , by (??) with c re- 
placed by cii. Note that only shear modes were con- 
sidered in |l6[ , an approximation which may hold for 
y <^ y* but misses the physics of the glass. Indeed 
only the compression modes are responsible for the glass 
(and lead to unbounded displacements for d > 2) since 
both displacements and force have to be considered along 
y. In d = 3 tilt modes would also be relevant for flux 
lines, and transverse shear modes for a solid. One finds 
y* = ciil{'nv)^J{cl^ - c2g)/(c2^ + cig). 

The predictions of glassy structure, topological order, 
channel motion and transverse critical force can be tested 
in experiments. For vortex lattices J}! can be measured 
in presence of a longitudinal current. Magnetic noise ex- 
periments and NMR probe vt -\- u{x,t) and the phase at 
the washboard frequency should contain a static compo- 
nent with slow and anisotropic decay. Other experimen- 
tal systems such as colloids, magnetic bubbles and dou- 
ble or triple incommensurate CDW should exhibit similar 
behavior. A transverse critical force may explain recent 
Hall voltage experiments in 2D Wigner crystals ||2^ . The 
predictions of channel motion and transverse critical cur- 
rent can be directly tested in numerical simulations [ p5| . 

(||) shows the importance of the relative orientation of 
the lattice and the applied force. It has been arg ued 
that to minimize power dissipation the lattice aligns with 
the force. This process may be slow and other orienta- 
tions can be studied by applying a transverse field. We 
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expect commensuration effects with devil's staircase type 
structure in the response to additional force. At higher 
commensuration vectors the channel structure may be- 
come unstable due to stronger effect of u'^. Conversely, 
the larger the static part, the more stable the glass with 
fewer topological defects. The size of the plastic flow 
regime should thus depend on the lattice orientation. It 
is then possible that in d = 3 for weak disorder the glass 
remains topologically ordered at all v and that the inter- 
mediate plastic regime disappears at low T. One would 
thus go smoothly from the moving to the static Bragg 
glass. At large v channels are nearly straight and out of 
equilibrium dislocations are thus suppressed. This allows 
to understand why in the plastic regime disappears 
when J is slowly decreased. 

Previous descriptions of moving systems, such as man- 
ifolds driven in periodic |Q or disordered potentials, fo- 
cused on the generation under RG of dissipative Kardar- 
Parisi-Zhang (Vu)^ non linearities. They do occur, due 
to lattice cutoff, in driven random sine-Gordon models 
|p7| , pT| . They are not important here because of the sta- 
tistical symmetry Uy — > —Uy in and the fact that 
dynamic modes are massive. Thus, because of peri- 
odicity this problem belongs to a new universality class. 
KPZ terms may play a role for incommensurate motion. 

In conclusion, we studied a lattice moving in a random 
potential. Static disorder dominates motion along sym- 
metry directions and the moving system is a glass with 
a large amount of topological order. It is continuously 
related to the static Bragg glass and although the decay 
of translational order is slow it has genuine glass proper- 
ties different from a usual solid. We predict experimental 
signatures such as elastic channels and transverse critical 
current. 
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